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Abstract 



^^ ' Total degree reverse lexicographic order is currently generally regarded as most 

C^ , often fastest for computing Grobner bases. This article describes an alternate less 

^ [ mysterious algorithm for computing this order using exponent subtotals and de- 

scribes why it should be very nearly the same speed the traditional algorithm, all 
other things being equal. However, experimental evidence suggests that subtotal 
^ I order is actually slightly faster for the Mathematica® Grobner basis implementa- 

lO ' tion more often than not. This is probably because the weight vectors associated 

with the natural subtotal weight matrix and with the usual total degree reverse 
lexicographic weight matrix are different, and Mathematica also uses those the 
^^ ' corresponding weight vectors to help select successive S polynomials and divisor 

^^ ■ polynomials: Those selection heuristics appear to work slightly better more often 

with subtotal weight vectors. 

However, the most important advantage of exponent subtotals is pedagogical. 
It is easier to understand than the total degree reverse lexicographic algorithm, 
p\ , and it is more evident why the resulting order is often the fastest known order for 

j^ I computing Grobner bases. 

Keywords: Term order, Total degree reverse lexicographic, tdeg, grevlex, Grobner basis 

1 Introduction 

Total degree reverse lexicographic order (degRevLex) is currently generally regarded as 
most often fastest for computing Grobner bases. HI This order is usually determined by 
Algorithm 1. However, as indicated in the comments therein, this algorithm tends to be 
mystifying and therefore difficult to recall: 
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^The adjective "graded" is sometimes used instead of "total degree" 



Algorithm 1 degRevLex order of two exponent vectors 



Given: Nonnegative integer exponent vectors a = [ai,a2, ■ ■ ■ , an], (3 = [/3i, /32, . . . , /^n]. 
Returns: One of "-<", "=", or 'V according to tlie degRevLex order between power 



products Zi^Z2'^ ■ ■ ■ z^" and z^^z^ ■ ■ ■ z^" witli indeterminate order Zi y Z2 y ■ ■ ■ >- 



Zr, 



a -^ ai; 
6^ A; 

for A; ^ 2 to n do 

a -i^ a + ak] 

b ^ b + (3k; end for; 
if a < 6, then return "-<"; 
if a > 6, then return 'V; 
for /c ■(— ??, to 1 by — 1 do 

if ak > f^k, then return "-<"; 

if ak < j3ki then return 'V; end for; 
return "="; 



/* . 



So far this all makes sense. */ 



/* : Whoa! Why is ">" matched with "^"? 
/* : These must be typographical errors!" 
/* : At least this step makes sense! */ 



This is the order that Buchberger in his Ph.D. dissertation and the order that Grobner 
always used when discussing multivariate polynomials (Buchberger, personal communi- 
cation) o 

The next fastest widely discussed order is total degree lexicographic order, and it too 
orders primarily by total degree. Therefore clearly total degree is very important for 
speed, and that makes sense because if we are iteratively annihilating terms with the 
largest total degrees, then the degree of each variable can't increase beyond that total 
degree. 

Consequently, it seems plausible that it would be more consistent to break total- 
degree ties with the sum of the degrees of all but the least main variable, then break 
those ties with the sum of the degrees of all but the two least main variables, and so on, 
as described in Algorithm 2. 



^Trinks [8] made the extraordinarily useful contribution of introducing the idea of admissible orderings 
and the particularly useful alternative example of lexicographic ordering. 



Algorithm 2 subtotal order of two exponent vectors 



Given: Nonnegative integer exponent vectors a = [ai,a2, ■ ■ ■ , an] , /3 = [/?i, /32, • • • , /?«]• 
Returns: One of "-<", "=", or ">-" according to the subtotal order between power prod- 
ucts z"^Z2^ ■ ■ ■ 2"" and Zi Z2 ■ ■ ■ z^" with indeterminate order Zi y Z2 >- ■ ■ ■ >~ Zn- 

Ai ^ tti; 
Bi ^ /?i; 
for A; ^ 2 to n do 

Ak ^ Ak-i + ak] 

Bk ^ Bk-i + 13k] end for; 
for A; ^ 77, to 1 by —1 do 

if Ak > Pk then return ">-"; /* : I have a better feeling about this! */ 

if Ak < Pk then return "-<"; end for; 
return "="; 

Algorithms 1 and 2 both do 2n additions followed by up to n comparisons, with the 
same looping costs. 

Both algorithms assume that the variables have been extracted from the power prod- 
ucts, which requires that all power products have the same number of exponents even if 
some of these exponents are 0. If instead the power products have only non-zero expo- 
nents, and therefore also contain variables to indicate the base for each exponent, then 
the algorithms are slightly different but both of them still have the same complexity as 
each other. 

Section [2] discusses weight matrices, weight vectors, and their implications for subtotal 
versus degRevLex order. Section [3] describes the experimental procedures for comparing 
these two orders and the results of those comparisons, with conclusions in Section |H 

2 Weight matrices and weight vectors 

When I first thought of subtotal order, 1 wondered if subtotal order would be faster 
than degRevLex order. Therefore 1 wanted a quick way to compare them experimentally 
without having to implement my own Grobner basis algorithm or learn the intricacies of 
an existing one to modify it. Fortunately some computer algebra systems permit users to 
specify their own ordering merely by providing a non-singular real square weight matrix 
W: To compare exponent vectors a and f3 for ordering, we lexicographically compare 
corresponding weight vectors 



w (a) ^— [wi {a) , W2 (ct) , . . 


.] ^ W-oi^ 


w{/3) ^ [w,{f3),W2{l3),.. 


. ^ W-0^ 



The weight matrix corresponding to algorithm 2 is 



W. 



sub 



/I 1 . 


• ^ ^\ 


1 1 . 


. 1 


1 1 . 


. 


\l . 


.00/ 



Proposition 1. The ordering specified by Wsnh is an admissible ordering. 

Proof. The first element of each column in Wsuh is positive. 

Also, matrix Wsuh is non-singular because the columns, hence the variables, can be 
permuted into an upper triangular matrix having 1 for every diagonal element, and the 
determinant of a square upper triangular matrix is the product of its diagonal elements, 
which consequently is non-zero, making the permuted Wsuh, hence also Wsuh non-singular. 
These are sufficient conditions for an admissible ordering. D 

Every admissible term order can be associated with a weight matrix, and the one that 
is usually associated with algorithm 1 for degRevLex ordering is 



W^ 



degRevLex 



/I 1 






1 1 \ 

-1 
-1 



(2) 



\ -1 ... y 

There are at least two ways to use weight matrices in an implementation: 

1. Do a matrix- vector multiplication for a power product every time we want to or- 
der it relative to another power product. To save some time, we could interleave 
the matrix-vector multiplication with the comparison of successive weight vector 
components to quit as soon as a difference is detected. 

2. Do matrix-vector multiplications only for the given input polynomials, store the 
weight vectors in parallel with the corresponding exponent vectors, then 

(a) whenever two power products are multiplied, compute the weight vector of 
their product as the sum of the two weight vectors, and 

(b) whenever two power products are divided, compute the weight vector of their 
quotient as the difference of the two weight vectors. 

At the expense of some additional storage space, the second choice is clearly much faster, 
because the number of power product comparisons during computation of a Grobner 
basis is usually far more than the number of power products in the given polynomials. 

Mathematica accepts weight matrices, and it uses the second of these two methods 
- even for the implemented built-in orderings, for which it uses essentially Algorithm 1 



rather than WdegRevLex to initiahze the weight vectors in the case of degRevLex ordering. 
Thus we can expect built-in degRevLex order for that implementation to be slightly 
faster than providing WdegRevLex, but usually not dramatically so. 

Nonetheless, to make the comparison fair - more as if comparing built-in degRevLex 
ordering to built-in subtotal ordering with initialization via Algorithm 2-1 started timing 
a few examples done with both PVgub and VFdegRevLex- I soon noticed that both matrices 
always gave the same Grobner basis. This could have been a coincidence. However: 

Proposition 2. Subtotal ordering is equivalent to degRevLex. 

Proof. The first rows of Wguh and WdegRevLex are identical. The second row of W^sub given 
by formula ([1]) is the second row of WdegRevLex given by formula ([2]) plus the first row 
of M/degRevLex- The third row of W^sub is the first row of IVdegRevLex plus the second and 
third rows of WdegRevLex, and so on. Thus IVdogRovLox can be transformed into Wguh by 
a sequence of adding positive multiples of rows to rows below them. This is a sufficient 
condition for orderings specified by two weight matrices to be equivalent. D 

Remark. After noticing the identical Grobner bases I discovered that Chee-Keng Yap [lOj 
describes the computation of degRevLex ordering by Algorithm 1, but lists Wgnh rather 
than the usual WdegRevLex as a weight matrix for degRevLex ordering. He didn't explain 
why, but I suspect that he had the idea of subtotals too. 

With weight matrices WdegRevLex and Wsub inducing the same ordering, we should 
expect their speed ratios to be close to LO. 

Section [3] indicates that this is often so, but not always, then explains why there are 
exceptions in the case of Mathematica. 

3 Experimental procedures and results 

The discussion in Sections [1] and [2] suggests that subtotal order should be very nearly 
the same speed as degRevLex order - regardless of whether they are both determined 
directly from exponent vectors or both determined via weight matrices and perhaps also 
weight vectors. This section describes experimental procedures and results that test this 
hypothesis. 

3.1 Experimental procedures 

A Mathematica function invocation of the form 

TimeConstrained [Timing [ 

GroebnerBasis [{polynomiali, polynomial2, ■ ■ ■} , {variablei, variable2, • • •} , 
Sort — )■ True, MonomialOrder — )■ DegreeReverseLexicographic]; ], 

maximums econds] 

computes a degRevLex Grobner basis of the polynomials after heuristically reordering 
the variables for speed, then displays only the computing time - or displays $ Aborted if 
maximums econds is exceeded. 



The rewrite rule 

SubtotalWeightMatrix [n_] := Table [Table [if [i < n - j + I, 1,0], j, n] , i, n] 

defines a function that returns an n by n weight matrix for subtotal ordering. 
The rewrite rule 

DegRevLexWeightMatrix [n_] := 

Table [Table [If [i == 1, 1, If [i == n + 2 - j, -1,0]], {j, n}], {i, n}] 

defines a function that returns an n by n weight matrix for subtotal ordering. 
Therefore a Mathematica function invocation of the form 

TimeConstrained [Timing [ 

GroebnerBasis [{poly nomiali, poly nomial2, • • •} , {variablei, . . . , variahltn} , 
Sort — )■ True, MonomialOrder — ^ SubtotalWeightMatrix [n]]; ], 

maximums econds] 

computes a subtotal Grobner basis of the polynomials after heuristically reordering the 
variables for speed, then displays only the computing time or $Aborted. A similar 
function invocation using DegRevLexWeightMatrix displays the computation time of 
a degRevLex basis using a weight matrix. 

I didn't want to address the issue of inexact computation at this time, because it 
is handled quite differently by different systemso Therefore I avoided test cases that 
contain Floats, unless they were obviously representations of rational numbers having 
small magnitude denominators, in which case I rationalized those Floats. 

I wanted to detect any difference in the relative speeds of term-order comparisons, 
and coefficient arithmetic tends to be a larger portion of the total computing time when 
Grobner bases are computed over the rational numbers rather than over the integers 
modulo a prime whose square fits in one computer word. Therefore, I did all problems 
in the coefficient domain ^32003, because the prime 32003 maps only about 0.003% of all 
non-zero integer coefficients to 0, but 32003 is small enough so that its square fits within 
one 32-bit computer word. 

Some of the original test cases and the ones obtained by rationalizing simple fioating- 
point coefficients contained rational coefficients that weren't integers. Such polynomials 
were multiplied by the least common multiple of their coefficient denominators so that 
computing a Grobner basis in Z32003 was straightforward. 

My patience for typing examples, checking for typographical errors, and waiting for 
results is limited. Also, I wanted to avoid the extra space of listing previously unpublished 
examples in two-dimensional .pdf format, forcing others to do lengthy error-prone typing 
to try all of these examples on some other Grobner basis implementation. Therefore 
I searched the Internet for medium-sized examples - preferably available in text that I 
could copy, paste and quickly edit to replace semicolons with commas, etc; and I used 120 



■^Many implementations make no special effort for Floats, making the results disastrously sensitive 
to differences in floating point arithmetic and differences in the order of operations. 

6 



seconds as the time limit. Many such examples are available at [9]. Others are available 
at [S] and fT] . These sites also give original references for the examples. This article lists 
input polynomials for a few additional examples that I think have not previously been 
published. 

Although the Sort — > True parameter should make the results rather insensitive to 
the order in {variablei, . . . , variablen}, I also entered the list of variables in the order 
specified in the sources or that I could infer, in case it mattered. Some of the original 
problems involved eliminating some of the variables or treating some as parameters in 
the coefficient domain. However, for uniformity I simply computed the Grobner basis 
with respect to all of the variables. 

My objective was to compare the speed of subtotal versus degRevLex ordering al- 
gorithms. To make the comparison fair, I used a weight matrix for degRevLex too. 
However, I also used the built-in degRevLex option to estimate how much improvement 
to expect if the subtotal ordering algorithm was built in. 

The computer has a 1.60GHz Intel Core 2 Duo U9600 CPU with 3 gigabytes of RAM. 

The Windows Vista operating system appears to have a timer resolution of only about 
0.015 seconds. Therefore if a time was less than 1 second, then I issued the command 

Timing [Do [ 

GroebnerBasis [{polynomiali, polynomial^, . . .} , {variablei, variablei, . . .} , 

Sort —7- True, MonomialOrder — ?■ . . .], {?ti}]; ], 

with the number of repetitions m sufficient to make the time exceed 1 second, then 
divided by m to obtain the time for one repetition. Despite this, even with no network 
connection and no voluntary programs launched other than Matheinatica, times tend 
to vary upon repetition within an interval of about ±3% of their mean. Mathematica 
uses reference counts rather than garbage collection, so this variation is probably caused 
instead by the multiple core architecture and the irregular competing activity of the many 
resident programs that lurk in most Windows installations. 

Relative speed is most important for problems that require extensive time. Therefore 
of all the problems that I tried, I included the ones that took the most time for built- 
in degRevLex order without exceeding the time limit for any ordering algorithm - as 
many examples as fit in a one-page table. The examples that I thus excluded didn't 
have noticeably different overall behavior regarding the relative speed of subtotal versus 
degRevLex ordering algorithms. 

3.2 The test examples 

Here are the included examples that I believe aren't already publicly published in some 
form: 



1. Lichtblau 10 

f^zb + x^ya, 

tx^yz — ab'^cde, 

xy^z'^d + zc^e^, 

tx'^y'^z'^ + ab'^c^e^, 

with variable order {t, x, y, z, a, b, c, d, e}. 

2. Lichtblau 2: 

a2 + 62 + 2c2 + 2d^ + 3/2 + 3g^ - h 



70 ab + 140 ac + 140 bd + 28 cd + 252 cf + 252 rf^ + 18/^ - 105 u, 

28 be + 28 c2 + 28 ad + 28 rf^ + 42 af + 12 df + 2Ap + A2bg+ 12 eg + 2Ag^- 35 w, 

36 ed + 30 bf + 24 c/ + 30ag + 24 rf^ + 16/^ - 35 x, 

8df + 2f + 8eg + 2g'-7y, 

100fg-77z, 

with variable order {a, b, c, d, /, (7, /i, u, v, x, y, z}. 

3. Lichtblau 3: 

-3375 uv + 3291 u\ + 750 uv^ - 732 ^/^^^ ^ 2225 uv^ - 2209 ^/^^jS + q^^ 

P + Qy, 
P + Qz, 

r (1 + 24 Mt; + 24 u^v - 24 uv^) , 

where 

P = -400 M+350 u^-4800 Mt;+4800 M^t;+7425 uv^-7359 u^v^-222'b Mt;^+2209 ?/V, 
Q = (\^2^uv^2^u^v-2^uv^), 

with variable order {x, y, 2, m, v, r}. 

4. A geometry problem of Michael Trott: 

-xx + 0:2 + ?/i - 3a;^?/i + 2x\yx - 2xxy\ + 2xiy\ -yi\ 3x\y2 - 2x\y2, 

x\-X2-y\^ 3x2?/i - 2x\yx + ?/2 - "ix^y-i + 22:^^2 + 2xi?/| - 2x2y\, 

-1 + 2x1 - 2x? + xf + 2yi + y\, 

-1 + 2x2 - 2x3 + x^ + 22/2 + ?/2, 

1 + (xi -X2)2; + (2/1 -?/2)^^ 

with variable order {xi, 3^2) l/i) ^21 -2}- 

5. Mathematica Help- GrobnerBasis- Options- Sort: 

3x'^ + 5xyz^ — 10 y'^z — 6xz + y^ + w, 

—2x'^z + 3x^y'^ + y^ — 12xz — 8xz^ + 3y'^z — 11 wxy"^, 

10 x'^w — 7yzw'^ — 2x2;'*u' + Ax'^y + 3xy'^ — 6yz^ — w + 2, 

w^ — wx'^y + xyz^ _ 2wxz'^ — 3w — 2xy'^ — 3, 

with variable order {w, x, y, z}. 



^The problems provided by Daniel Lichtblau are from a mix of literature, user questions and bug 
reports, with unrecorded individual provenance. 
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6. Variation on a theme of Giovini et. al: One of my test files is so similar to Giovini 
3.7 [4j that I must have entered it twice, but once with typographical errors. These 
seemingly minor changes make the example that is relatively nearly fastest for 
subtotal ordering algorithm become an example that is relatively slowest for that 
algorithm. The variation is 

^33^23 _ y82^^ 
^45 _ ^13^21^^ 

x'^^c — y^^z^"^, 

2;22 _ ^332-12^^ 
^5^17^22g _ ^^ 

xyzt — 1, 
with variable order {t, b, c, e, d, a, z, x, y}. 

All of the other examples are already published elsewhere, as cited in Table 1. 

3.3 Test results 

The last column of Table 1 displays the most important results - the ratios of the com- 
puting time for subtotal ordering versus degRevLex ordering, both using a weight matrix. 

Notice that the rows are ordered by non-decreasing values in this last column in an 
attempt to discern correlations with the number of variables and/or total degrees of the 
input polynomials. These total degrees are also listed in the Table with, for example, 
6^-5 meaning 3 polynomials each having total degree 6 and 1 polynomial having total 
degree 5. However, the number of examples and/or their variability doesn't appear to be 
large enough to reveal obvious correlations. This might be partly because the standard 
deviation of the number of variables is only 2.8 for a median and mean of about 8. 

Most of the speed ratios are rather close to 1.0, as expected. However, there are 
some notable outliers at at the top and bottom of the table. A probable explanation for 
these outliers is that in Mathematica the weight vectors are also used to select the next 
S-polynomial or reducer polynomial. If so, then it is an indication that good selection of 
a next S-polynomial and a next reducer polynomial should instead if possible be based 
on the ordering that the weight vector induces, so as to be invariant to a property that 
doesn't correlate completely with ordering. 

Although re-executing an example caused time variations for individual examples 
within an interval of about ±3%, there are enough examples so that the summary statis- 
tics are more tightly repeatable. 

The median speed ratio of 0.98 and mean of 0.92 suggest that the subtotal algorithm 
would be about this much faster than the degRevLex algorithm if both were built into 
Mathematica. The standard deviation of 0.24 weakens that conclusion, but most of the 
standard deviation comes from the outliers at the top of the table where the subtotal 
algorithm was significantly faster. 

The speeds of the two methods are close enough so that quite possibly the degRevLex 
algorithm could be faster than the subtotal algorithm if both were built into another 
Grobner basis implementation that used weight vectors in a different way to select S 



polynomials and reducers. However, it seems highly likely that the speeds would be 
quite close if weight vectors weren't used for selection strategy. 
The penultimate column of Table 1 is next most interesting: 

• The summary statistics for that column weakly support the conclusion that built-in 
degRevLex is faster than the degRevLex weight matrix by only 2% for the median 
or 8% for the mean. This can only be attributable to computing the initial weight 
vectors from the weight matrices, and these small percentages indicate that the 
initialization of weight vectors is usually only a small portion of the total time. 

• However, there are a surprising number of instances where the weight matrix is 
slightly faster - up to 7%. I can think of no implementation reason why this 
should be. If there is no such reason, then it can be taken as an indication of the 
repeatability deviations in individual time ratios - up to 7% rather than the 3% 
that I estimated from informal experiments. 



4 Conclusions 

The main result is that subtotal ordering is an alternate way to view degRevLex ordering 
that more clearly explains its good behavior. 

A secondary conclusion is that the the algorithm for computing subtotal order is very 
nearly the same speed as that for computing degRevLex order - at least if the induced 
ordering rather than the weight vectors is used to select the next S-polynomial and next 
reducer. 

The relative speeds of the subtotal and degRevLex ordering algorithms are close 
enough so that it probably isn't worth replacing the degRevLex algorithm with the 
subtotal algorithm in existing implementations. However, it is well worth considering use 
of the subtotal algorithm instead of the degRevLex algorithm in new implementations. 
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Table 1: Coefficient in Z,32oo3- Seconds & time ratios for subtotal vs degRevLex order 


search term & citations 


# 

vars 


input tot. degrees 


grevlex 

builtin 

sec. 


grevlex 

/ grevlex \ 
V matrix / 


subtotal 

^ grevlex \ 
V matrix / 


CohnS UJ 


4 


63-5 


7.24 


0.08 


0.08 


filter design [6] 


9 


5, 42.3, 2^ 


10.1 


0.98 


0.22 


benchmark il [IJ 


10 


310 


2.78 


0.30 


0.31 


Assur44 |lj 


8 


3^-25 


6.30 


0.40 


0.38 


Giovini 3.7 [4J 


9 


83-453-44-4 


36.7 


0.95 


0.65 


benchmark Dl yj 


12 


32-29-1 


0.71 


0.59 


0.73 


des22_24 P E] 


10 


28.12 


0.75 


0.70 


0.73 


Lichtblau 2 


9 


11-10-62 


0.44 


0.72 


0.90 


Gonnet et. al. [5j 


17 


219 


6.01 


0.74 


0.95 


cdpmS [Ij 


5 


35 


4.60 


0.95 


0.95 


reimerS [9j 


5 


6-5-4-3-2 


1.59 


0.99 


0.96 


Kotsireas4body [6j 


6 


53.23 


2.92 


1.05 


0.97 


Lichtblau 3 


12 


2« 


0.50 


0.88 


0.97 


cyclic6 [1, 9J 


6 


6-5-4-3-2,-l 


0.62 


1.01 


0.97 


Giovini 3.1 [4j 


7 


42. 310. 2 


0.57 


1.17 


0.97 


eco8 [Ij 


8 


33-2-1 


1.18 


0.95 


0.98 


redeco7 [9j 


8 


26.12 


0.67 


0.96 


0.98 


extcycS [Ij 


6 


52-4-3-2-1 


1.15 


1.00 


0.98 


f744 [IJ 


12 


32.22.12 


5.13 


0.93 


0.98 


Lichtblau 1 


6 


53 


5.91 


1.02 


0.99 


virasoro [9j 


8 


2« 


23.8 


0.99 


1.00 


Trott geometry 


5 


43.3 


12.4 


0.99 


1.00 


Kotsireaus4bodySymmetric f6\ 


7 


52-32-22 


36.2 


1.00 


1.01 


Katsura7 [9J 


7 


2^.1 


0.78 


0.99 


1.01 


redcyc6 [9j 


6 


11-5-4-3-2-1 


0.43 


1.00 


1.01 


Harrier RK2 [?] 


13 


44.32.2.14 


33.5 


0.95 


1.03 


Mathematica help 


4 


7-6-5-4 


5.98 


1.01 


1.04 


rpbl24 [9J 


9 


32-26.1 


1.87 


1.07 


1.04 


Tran, rational implicitization [7J 


5 


62-5 


5.13 


0.96 


1.05 


kinema [Ij 


9 


26.13 


2.07 


1.03 


1.07 


KotsireausSbody [6, 9j 


6 


53.23 


3.18 


1.03 


1.10 


rpbl [9J 


6 


3^-2 


0.94 


1.07 


1.13 


variation on Giovini 3.7 


9 


83.46.453.4 


19.8 


0.94 


1.39 


Statistics! 












Median 


8.0 




2.85 


0.98 


0.98 


Mean 


8.1 




7.34 


0.92 


0.92 


Standard Deviation 


2.8 




10.7 


0.19 


0.24 


Count < 1.0000 








22 


20 


Count > 1.0000 
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